R Notebook

library(tidyverse)
library(sf)
library(rmapshaper)
library(here)
library(leaflet)
library(htmltools)

scotland <- read_sf(here("data_raw/SG_NHS_HealthBoards_2019/SG_NHS_HealthBoards_2019.shp"))

class(scotland)

head(scotland, 14)

Simple feature collection with 14 features and 4 fields
Geometry type: MULTIPOLYGON
Dimension:     XY
Bounding box:  xmin: 5512.998 ymin: 530250.8 xmax: 470332 ymax: 1220302
Projected CRS: OSGB 1936 / British National Grid

scotland_geometry <- st_geometry(scotland)

scotland_geometry

scotland_simple <- ms_simplify(scotland)

plot(scotland_simple)

obesity <- read_csv("../data_clean/obesity.csv") %>%
  filter(sex == "All") %>%
  filter(date_code == "2016-2019") %>%
  filter(obesity == "Obese")

scotland_simple %>%
  left_join(obesity, by = c("HBCode" = "feature_code")) %>%
ggplot() +
  geom_sf(aes(group = hb_name, fill = value)) +
  theme_minimal() +
  theme(panel.grid.major = element_blank(), axis.text = element_blank()) +
  ggtitle("Obesity in Scotland") +
  scale_fill_gradient(low = "green", high = "red")

smoking_status <- read_csv("../data_clean/smoking_status.csv") %>%
  filter(sex == "All") %>%
  filter(date_code == "2016-2019") %>%
  filter(smoking_status == "Current smoker")

scotland_simple %>%
  left_join(smoking_status, by = c("HBCode" = "feature_code")) %>%
ggplot() +
  geom_sf(aes(group = hb_name, fill = value)) +
  theme_minimal() +
  theme(panel.grid.major = element_blank(), axis.text = element_blank()) +
  ggtitle("Smoking Status in Scotland") +
  scale_fill_gradient(low = "green", high = "red")

e_cig <- read_csv("../data_clean/e_cig.csv") %>%
  filter(sex == "All") %>%
  filter(date_code == "2016-2019") %>%
  filter(e_cigarette_use == "Currently using")

scotland_simple %>%
  left_join(e_cig, by = c("HBCode" = "feature_code")) %>%
ggplot() +
  geom_sf(aes(group = hb_name, fill = value)) +
  theme_minimal() +
  theme(panel.grid.major = element_blank(), axis.text = element_blank()) +
  ggtitle("E-cigarette use in Scotland") +
  scale_fill_gradient(low = "green", high = "red")

alcohol %>%
  filter(sex == "All") %>%
  filter(hb_name == "Scotland") %>%
    filter(date_code == "2016-2019") %>%
  filter(alcohol_consumption == "Hazardous/Harmful drinker")

alcohol_scotland


  select(value) %>%
  pull()

obesity_scotland <- read_csv("../data_clean/obesity.csv") %>%
  filter(sex == "All") %>%
  filter(hb_name == "Scotland") %>%
    filter(date_code == "2016-2019") %>%
  filter(obesity == "Obese")

obesity_scotland

alcohol %>%
                    filter(sex == "All" | date_code == "2016-2019") %>%
                    filter(hb_name == "Scotland") %>%
                    filter(alcohol_consumption == "Hazardous/Harmful drinker") %>%
                    select(value) %>%
                    pull()

smoking_status_scotland <- read_csv("../data_clean/smoking_status.csv") %>%
  filter(sex == "All") %>%
  filter(hb_name == "Scotland") %>%
    filter(date_code == "2016-2019") %>%
  filter(smoking_status == "Current smoker")

smoking_status_scotland

e_cig_scotland <- read_csv("../data_clean/e_cig.csv") %>%
  filter(sex == "All") %>%
  filter(hb_name == "Scotland") %>%
    filter(date_code == "2016-2019") %>%
  filter(e_cigarette_use == "Currently using")

e_cig_scotland

unique (obesity$sex)

                selectInput("year_selection",
                            "Year Selection",
                            choices = unique(obesity$date_code))

alcohol %>%
                    filter(sex == "All" , date_code == "2016-2019") %>%
                    filter(hb_name == "Scotland") %>%
                    filter(alcohol_consumption == "Hazardous/Harmful drinker")